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ABSTRACT 

Brown dwarf formation and star formation efficiency are studied using a 
nested grid simulation that covers five orders of magnitude in spatial scale (10 4 - 
0.1 AU). Starting with a rotating magnetized compact cloud with a mass of 
0.22 M , we follow the cloud evolution until the end of main accretion phase. 
Outflow of ~ 5 km s -1 emerges ~ 100 yr before the protostar formation and does 
not disappear until the end of the calculation. The mass accretion rate declines 
from ~ lO _6 M yr" 1 to ~ 1O~ 8 -1O~ 12 M yr" 1 in a short time (~ 10 4 yr) after 
the protostar formation. This is because (1) a large fraction of mass is ejected 
from the host cloud by the protostellar outflow and (2) the gas escapes from 
the host cloud by the thermal pressure. At the end of the calculation, 74% 
(167Mj up ) of the total mass (225Mj up ) is outflowing from the protostar, in which 
34% (77Mj up ) of the total mass is ejected by the protostellar outflow with super- 
sonic velocity and 40% (90M Jup ) escapes with subsonic velocity. On the other 
hand, 20% (45Mj up ) is converted into the protostar and 6% (13Mj up ) remains 
as the circumstellar disk. Thus, the star formation efficiency is e = 0.2. The 
resultant protostellar mass is in the mass range of brown dwarfs. Our results 
indicate that brown dwarfs can be formed in compact cores in the same manner 
as hydrogen-burning stars, and the magnetic field and protostellar outflow are 
essential in determining the star formation efficiency and stellar mass. 
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— MHD — stars: formation — stars: low-mass, brown dwarfs 
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Introduction 



There are several scenarios for the formation of brown dwarfs (see, IWhitworth et al. 



20071 ). though many observations indicate that brown dwarfs are for med through the gravi- 



young brown dwar: 



2003h . Recently, 



s show the existence of circumstellar disks (e.g. 



Luhman et al.ll2005f) and of gas accret ion onto the young brown dwarfs (e.g.. lWhite fc Basri 



Pascucci et al. 


20031; 


(e.g., 


White & Basri 



Whelan et al. (boosh observ e d an optical jet around a young brown dwarf 



of p Oph 102. In addition, iPhan-Bao et al.l (120081 ) observed molecular outflows from the 
same object. To dat e, several outflows have been observed around young brown dwarfs 
( jWhelan et al.l l2009al ibl). Since outflow is typical in the star formation process, it is con- 



sidered that these flows are direct proof that brown dwarf formation occurs through the 
gravitational collapse of a molecular cloud core. On the other hand, observations indicate 
that the prestellar core mass function (CMF) resembles the initial mass function (IMF, e.g., 
Motte et al.lll998l ). suggesting that the effective reservoirs of mass required for the forma- 
tion o f individual stars in cluding brown dwarfs are already determined at the prestellar core 
stage (lAndre et al.l 120081 ) . The similarity of CMF and IMF implies that the prestellar cores 
are on a one-to-one correspndence with protostars, in which the star formation efficiency is 
expected to be e = M star /M core ~ 20-50%. 

Not all the mass of a prestellar core can convert into a star, be cause the outflow from a 
protostar can remo v e a si gnificant mass from the natal cloud core (IWhitworth et al.l 120071 ) . 
Matzner fc McKed ( 120001 ) has pointed out that since the gas accretion onto a protostar 
is impeded by the protostellar outflow, only a fraction of the mass can convert into the 
protostar, with the remainder being blown out. They modeled the prestellar outflow and 
analytically estimated a star formation efficiency of e ~ 25-70%. 

We can presume that brown dwarfs are born in more compact (or less massive) cloud 
cores, because they are less massive than stars. Taking a star formation efficiency of e = 30%, 
brown dwarfs with masses of 13-75 Mj up , where Mj up is the Jovian mass, can be formed in 
cores with masses of 43-250 Mw Cloud cores with such masses have been ob served in 



several surveys (e.g.. iMotte et al.ll200ll ; ISandell fc Kneel 1200 ll ; lOnishi et al.ll2002t ). In this 



paper, we study the evolution of a compact molecular cloud core and the formation of a 
brown dwarf. 
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2. Model and Numerical Method 

To study the evolution of star-forming cores in a large dynamic range of spatial scale, 
a three-dimensional nested grid method is used, in which the equations of resistive magne- 
tohydro dynamics wi th a barotropic equation of state are solved (see Eq. [l]-[3] and [5] of 



Machida et all 120071 ) 



In the collapsing cloud core, we assume protostar (or proto-brown dwarf) formation 
occurs when the number density exceeds n > 10 14 cm -3 at the cloud center. To model 
the protostar, we adopt a sink around the center of the computational domain. In the 
region r < r sin k = 0.5 AU, gas having a number density of n > 10 14 cm -3 is removed from 
the computational domain and added to the protostar as a gravity in each timestep. In 
addition, inside the sink, the magnetic flux is removed by the Ohmic dissipation (The detailed 
description will appear in a subsequent paper). 

As the initial state we take a spherical cloud with a Bonnor-Ebert (BE) density profile 
that extends up to twice the critical BE radius (r = 2 R c ) . Outside the sphere (r > 2 R c ) , a 
uniform density is adopted. For the BE density profile, we adopt a central density of n c = 
5 x 10 7 cm -3 and an isothermal temperature of T = 10 K. For these parameters, the critical 
BE radius is R c = 650 AU. Thus, the radius of the initial sphere is r = 2 R c = 1300 AU. To 
promote the contraction, we increase the density by a factor of 1.68. The density contrast 
between the cloud center and ambient medium is about 80. The mass within r < 2 R c 
(r < R c ) is 0.22 M (0.11 M Q ) corresponding to 225M Jup (112M Jup ). In this paper, we call 
the sphere with r < 2 R c the host cloud. The cloud rotates rigidly with Q = 2.5 x 10~ 12 s _1 
around the z-axis in the region of r < 2 R c , while a uniform magnetic field (B Q = 212 /xG) 
parallel to the z-axis (or rotation axis) is adopted over the whole computational domain. 
In the region of r < 2 R c (r < R c ), the ratios of thermal «o, rotational fio, and magnetic 
7o to the gravitational energy are a = 0.42 (0.5), (3 = 0.04 (0.016), and 7 = 0.13 (0.04), 
respectively. The gravitational force is ignored outside the host cloud (r > 2 R c ) to mimic a 
stationary interstellar medium. 



sec 



To calculate oyer a la rge spatial scale, the nested grid method is adopted (for details 



Machida et al. ll2005a Jbl). Each level of a rectangular grid has the same number of cells 
(64 x 64 x 32). The calculation is first performed with five grid levels (I = 1-5). The box 
size of the coarsest grid I — 1 is chosen to be 2 5 R C . Thus, a grid of I = 1 has a box size 
of ~ 2 x 10 4 AU. A new finer grid is generated before the Jeans condition is violated. The 
maximum level of grids is restricted to / max ^ 10. The I = 10 grid has a box size of 40 AU 
and cell width of 0.63 AU. 
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Results 



Since we adopted an unstable core with a centrally peaked density profile as the initial 
state, the central part of the cloud collapses first and the gas density increases with time. 
In the collapsing cloud core, just after the cen tral number density reaches n r ~ 10 11 cm ~ 3 , 



the first adiabatic core (hereafter the first core, lLarsonlll969l ; iMasunaga fc Inutsuka Il2000l ) is 
formed with shock 8.35 x 10 3 yr after the cloud collapse begins. The first core has a disk-like 
shape at its formation, with a size of 18 AU in the cylindrical radial direction and 3.2 AU 
in the vertical direction with a mass of 1.1 x 10~ 2 M©. A low-velocity outflow (hereafter, 
outflow) of ~ 2 kins -1 appears around the center of the collapsing cloud 121 yr after the first 
core formatio n. This kind of outflow in a collapsing cloud has been reported in many other 



studies (e.g.. iTomisakal 120021 ; iBanerjee &: Pudritzl 120061 ; iHennebelle &: Fromangi 120081 ) . The 



central density exceeds n c > 10 cm and the protostar (hereafter, the proto-brown dwarf) 
is formed 242 yr after the first core formation. Thus, the outflow begins to be driven ~ 100 yr 
before the proto-brown dwarf formation. We calculated the cloud evolution 2.3 x 10 4 yr after 
the proto-brown dwarf formation, in which the first core increases in size keeping a disk-like 
structure and smoothly becomes the circumstellar disk (hereafter the circum-brown dwarf 
disk) with time. The outflow continues to be driven from the circum-brown dwarf disk and 
extends up to ~ 5 x 10 3 AU with a maximum speed of 7.2 km s -1 . From the start of the 
cloud collapse we calculated the cloud evolution for 2.9 x 10 4 yr, corresponding to 11.7 tg t0 , 
where tf^o is the freefall timescale at the center of the initial host cloud. 

Figure [1] shows the mass accretion rate for the proto-brown dwarf, and the mass of the 
proto-brown dwarf, outflowing gas, and circum-brown dwarf disk against time (i c ) after the 
proto-brown dwarf formation. The outflowing gas is defined as gas with velocity v r > c s 
for the whole computational domain, while the disk is defined as gas with v r < c s and 
n c > 5 x 10 12 cm -3 inside the host cloud. Figure [1] shows that the mass accretion rate 
remains almost constant at M ~ lO _6 A / / yr _1 for t c < 1.5 x 10 4 yr, while it suddenly drops 
and has a very small value of M ~ 1O~ 8 -1O~ 12 M yr -1 for t c > 1.5 x 10 4 yr. Note that it 
slightly decreases with time even for t c > 1.5 x 10 4 yr. Reflecting the accretion rate, the mass 
of the proto-brown dwarf continues to increase for t c < 1.5 x 10 4 yr, while it remains constant 
for t c > 1.5 x 10 4 yr. At the end of the calculation, the mass of the proto-brown dwarf is 
M ps = 45Mj up . Since the mass accretion rate is very small for t c > 1.5 x 10 4 yr, this object 
cannot acquire enough mass to become a hydrogen-burning star and hence evolves into a 
brown dwarf. On the other hand, the disk mass continues to increase for t c < 1.9 x 10 4 yr 
and saturates at Mdi S k ~ 10Mj up . Finally, the disk mass reaches Mdi S k = 13Mj up . Even if all 
the mass of the disk falls into the central object, the central object is within the brown-dwarf 
mass range (M ps + M disk = 58M Jup ). 
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Figure [T] shows that the mass of outflowing gas increases over time and reaches M out = 
128Mj up . Although outflow weakens for t c > 2 x 10 4 yr, it continues to be driven from the 
disk until the end of the calculation. Figure [2] shows the configuration of the outflow. The 
iso-velocity of v z = c s is represented by the transparent red surface, inside which the gas 
is outflowing from the center of the cloud with supersonic velocity. The red surface at the 
center indicates an iso-density of n = 2 x 10 6 cm" 3 , which almost corresponds to the density 
at the border between the initial host cloud and the envelope (n amp = 1.1 x 10 6 cm~ 3 ). The 
sphere enclosed by the white dotted circle (r = 2r cri ) represents the initial host cloud. The 
figure shows that the outflow penetrates the host cloud to reach ~ 5 x 10 3 AU, about four 
times larger than the initial host cloud (1300 AU). 

Figure [2] shows that the outflow has a wide opening angle inside the host cloud, while 
it has good collimation outside the host cloud. An hourglass-like configuration of field lines 
is realized inside the host cloud, because the field lines converge toward the center as the 
cloud collapses. The gas flows along the hourglass-like lines inside the host cloud and the 
outflow has a wide opening angle. On the other hand, gas flows along the straight field lines 
and has good collimation outside the host cloud. The opening angle of the outflow inside 
the host cloud strongly influences the mass accretion rate and star formation efficiency 



flMatzner fc McKedl2000h 



The density and velocity distribution around and inside the host cloud are plotted in 
Figure [3J Figure [3]a indicates that the residual matter in the host cloud is mainly distributed 
in a region along the rotation axis (i.e., z-axis) and on the equatorial (z = 0) plane. The 
source of the matter near the rotation axis is the outflow, while the matter on the z = 
plane corresponds to the disk which is supported by the centrifugal force. At the end of the 
calculation, M res = 104M Jup of the residual mass remains inside the host cloud (r < 2R C ). 
Since the mass of the proto-brown dwarf is M ps = 45Mj up and the initial host cloud mass is 
M ini = 225M Jup , the mass ejected from the host cloud is M e j = M ini — M rcs — M ps = 76M Jup . 
In addition, since the outflowing mass is M out = 128Mj up , the mass swept by the outflow into 
interstellar space is M sw = M out — M e j = 52Mj up . Thus, in the host cloud, 34% (= M e j/M in i) 
of the total mass is ejected by the outflow, while 20% (= M ps /M ini ) is converted into the 
star (or proto-brown dwarf). 

In addition, at the end of the calculation, a large fraction (M esc = 88Afj up ) of the residual 
mass (M res = 104Mj up ) has a positive flow velocity (v r > 0) and escapes from the host cloud. 
Note that we defined the outflowing mass M out as the gas with v r > c s in Figure [fl while 
we defined the escaping mass M esc as the gas with v r > inside the host cloud. The rest 
of the mass (i.e., the gas with v r < 0) inside the host cloud is 16Mj up (= M res — M esc ), 
which is composed of two parts: the circum-brown dwarf disk and the accreting gas. As 
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denoted above, since the disk mass is Mdisk = 13Mj up , the mass of the accreting matter is 
= 3M Jup . Thus, the mass of the accreting matter is 3% (= M acc /M res ) of the residual 
mass [or 1% (= M acc /M ini ) of the mass of the initial host cloud]. As a result, it is considered 
that the gas accretion is nearly finished. However, since the gas supply into the disk does 
not completely halt, (weak) outflow continues until the end of the calculation. 

In Figure [3]a, we divided the host cloud into three zones, according to our results. The 
outflowing gas has a supersonic velocity (v r > c s ) in the outflowing zone, while it has a 
subsonic velocity (0 < v r < c s ) in the escaping zone. Note that almost all the matter in 
the escaping zone exceeds the escape velocity of the host cloud. A part of the matter in the 
escaping zone comes from outflow from the disk, while the matter near the border between 
the inflow and outflow (i.e., near the orange contour of v r = 0) escapes from the host cloud 
by the pressure gradient force. Since 34% of the total mass is already ejected from the host 
cloud, the gravity inside the host cloud weakens. As the initial state, we adopt a nearly 
equilibrium state, in which the gravity is balanced with the thermal pressure gradient force. 
Thus, owing to a decrease of mass (or gravity), the matter near the envelope escapes as the 
outflow continues. While the outflow remains inside the host cloud, the gas in the escaping 
zone in Figure [3]a has a negative radial velocity (v r < 0) and falls toward the proto-brown 
dwarf. After the outflow penetrates the host cloud and a large fraction of mass is ejected, 
the radial velocity gradually decreases and finally becomes positive (v r > 0), thereafter the 
gas escapes. 

Figure and [3]c are close-up views of Figure [3]a. As seen in Figure strong outflow 
appears in the region of r < 100 AU. In addition, a cavity wall appears around the border 
between the outflow and inflow. During the calculation, strong outflow continues intermit- 
tently. Figure[3]c shows the roots of the outflow located at ~ 10 AU, far from the proto-brown 
dwarf. 



4. Discussion 



Recently, iPhan-Bao et al.l (120081 ) observed bipolar molecular outflow from a proto-brown 
dwarf (p Oph 102) with ~ 60Mj up . The outflow extends up to ~ 1000 AU with a maximum 
velocity of 2.2 km s -1 and mass of 0.17Mj up . In addition, they estimated that the circum- 
brown dwarf disk has a size of ~ 80 AU and mass of 8.3Mj up . In our calculation, the outflow 
driven from the proto-brown dwarf with 45Mj up has a typical flow speed of ~ 2kms~ 1 
(Fig. [3f)) and mass of 128Mj up , and extends up to ~ 5 x 10 3 AU. The circum-brown dwarf 
disk in our calculation has a mass of M^sk = 13Mj up and size of ~ 150 AU. Thus, our results 
are quantitatively consistent with observations, except for the outflowing mass. 
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Phan-Bao et al.l (120081 ) and IWhelan et al.l (120051 ) estimated the mass accretion rate to 
be M ~ lO _9 M yr _1 , which is comparable to our result in the later accretion phase (t c > 
1.5 x 10 4 yr). However, their scenario for brown-dwarf formation is significantly different 
from our results. We consider that a part of a brown dwarf is formed in a (small, compact) 
cloud core in a similar process as low-mass star formation. However, in observational studies 
it seems to be considered that the brown dwarf formation is a scaled down version of that of 
hydrogen-burning stars, in which a very small rate of mass accretion (M ~ 10~ 9 -10~ 12 M & ) 
is assumed in the main accreton phase. We show that, even in the brown dwarf formation 
process, the mass accretion rate remains as high as M ~ 1O _6 M yr -1 for t c <^ 10 4 yr, and 
it rapidly drops to M ~ 1O" 9 -1O" 12 M yr" 1 for t c ^ 10 4 yr. 

In the theoretical framework of star formation, the mas s accretion rate can be d escribed 
as M = fcl/G, where / is a constant (e.g., / = 0.975 for IShulll977l . / = 46.9 for iHunter 
19771 ). Since gas clouds have temperatures of T ^ 10 K (c s ^ 0.2 kms -1 ), the accretion rate 
in the main accretion phase is M ^ 10 _6 M Q yr _1 which is the minimum value attained in 
the general star formation scenario. Our results show, in the main accretion phase (t c < 
1.5 x 10 yr), the accretion rate of ~ 2 x 10~ 5 — 1O~ 6 M yr -1 that corresponds to / ~ 1 — 10. 
Note that other effects such as cloud rotation, magnetic field, and turbulence only increase 
the mass accretion rate. Thus, the scenario in which the accretion rate of 1O~ 9 M yr -1 
lasts for ~ 10 7 yr is not theoretically expected. On the contrary, it is reasonable that an 
accretion rate of ~ 1O~ 6 M yr -1 lasts for ~ 10 4 yr for brown-dwarf formation. In this case, 
since the main accretion phase is very short, it is difficult to observe the accretion onto a 
proto-brown dwarf with M ~ lO~ 6 Af yr -1 . Thus, we may often observe low accretion rates 
of < 1O~ 9 M yr -1 , because a lower accretion phase is expected to last for a longer duration 
(3> 10 4 yr). In this study, the accretion rate drops at t c ^ 10 4 yr because we adopted a less 
massive core as the initial state. However, we expect that if a massive core is adopted as the 
initial state, a high accretion rate of ~ 10~ 6 M Q yr _1 lasts for a long duration of t c ~ 10 5 - 
10 6 yr, and a stellar mass star forms. The result in this letter favors the idea that the stellar 
mass is determined by the initial size of the host cloud and hence brown dwarf sized objects 
are formed in less massive clouds. 

We comment on the sudden drop of the mass accretion rate at t c ~ 10 4 yr, as seen 
in Figure [TJ The epoch of the sudden drop almost corresponds to the epoch at which the 
outflow escapes from the host cloud. Since the outflow with a significant mass is escaped 
from the host cloud by this epoch, the gravitational potential of the host cloud shallows and 
the accretion rate weakens, which causes the sudden drop of the accretion rate. We expect 
that the epoch of this sudden drop depends on the size of the host cloud. 



For the outflow model, iPhan-Bao et al.l (120081 ) proposed a jet-driven bow shock model 
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( IMasson &: Cherninl Il993l ). However, our results are different from the predictions of this 
model. Our calculations show that the outflow is directly driven from the circum-brown 
dwarf disk. No appearance of the high-speed jet in our calculation is due to the fact that 
we adopted the region of r < 0.5 AU as the sink, and thus, we did not calculate the region 
near the proto-brown dwarf. Adopting a central object with a mass of 40Mj up , the Kepler 
speed at ~ 0.5 AU is ~ lOk ms" 1 . Thus, f l ow of > lOkms -1 could not be resolved in 
our calculation. In contrast, iMachida et al.l (j2008bl ) calculated the structure in the near 
proximity of the protostar and showed that a high-speed jet of ~ 30kms _1 is driven from 
the Jovian-mass protostar. Thus, when we resolve the structure in the near proximity to the 
proto-bro wn dwarf, a high-speed jet is expected to be driven from the proto-brown dwarf. As 
shown in IMachida et al.l (j2008bl ). we expect that the high-speed jet hardly affects the mass 
accretion rate and star formation efficiency, because it has a well-collimated structure. Our 
calculations show that the wide-opening molecular outflow driven from the circum-brown 
dwa rf disk strongly af f ects th e mass accretion rate and star formation efficiency, as predicted 



in 



Matzner & McKeel (120001 ) . 
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Fig. 1. — Mass accretion rate (diamonds) , and mass of the proto-brown dwarf (red), circum- 
brown dwarf disk (green), and outflowing gas (blue) against the elapsed time after the proto- 
brown dwarf formation. The left axis indicates the mass accretion rate, while the right axes 
indicates mass in units of Jovian and solar masses. The two horizontal dashed lines indicate 
the lower and upper limits of brown dwarf mass. 



Fig. 2. — Three dimensionsal structure in the second largest grid (1 = 2). The structure of 
the outflow (v z > c s ) is shown by the transparent red surface. The structure of the high- 
density region (n = 2 x 10 6 cm -3 ; red-isodensity surface) and magnetic field lines (black and 
white streamlines) are plotted. The density contours (colors) on the x = plane and velocity 
vectors (arrows) on the y = plane are projected on each wall surface. The white dotted 
line and transparent black sphere show the initial host cloud of r = 2R C . 
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Fig. 3. — Density {color scale) and velocity (arrows) distribution on the y = plane at 
t = 12.71 t ff (~ 3.18 x 10 4 yr) for different grids of I = 3 (a), 6 (b), 9 (c). The orange 
contour in each panel corresponds to the boundary between the outflowing (t> r > 0) and 
inflowing (u r < 0) gas. The red contour in panel (a) indicates v z = c s , inside which the 
gas is outflowing with supersonic speed. The white dotted circle in panel (a) indicates the 
radius of the initial host cloud r = 2R C . The sphere inside r < 2R C (i.e., inside the white 
dotted circle) is divided into three regions: the accreting [y r < 0), escaping (0 < v r < c s ) 
and outflowing {v r > c s ) zones. 



